Alu RNA fold links splicing with signal recognition particle proteins.

Transcriptomic diversity in primates was considerably expanded by exonizations of intronic Alu elements. To better understand their cellular mechanisms we have used structure-based mutagenesis coupled with functional and proteomic assays to study the impact of successive primate mutations and their combinations on inclusion of a sense-oriented AluJ exon in the human F8 gene. We show that the splicing outcome was better predicted by consecutive RNA conformation changes than by computationally derived splicing regulatory motifs. We also demonstrate an involvement of SRP9/14 (signal recognition particle) heterodimer in splicing regulation of Alu-derived exons. Nucleotide substitutions that accumulated during primate evolution relaxed the conserved left-arm AluJ structure including helix H1 and reduced the capacity of SRP9/14 to stabilize the closed Alu conformation. RNA secondary structure-constrained mutations that promoted open Y-shaped conformations of the Alu made the Alu exon inclusion reliant on DHX9. Finally, we identified additional SRP9/14 sensitive Alu exons and predicted their functional roles in the cell. Together, these results provide unique insights into architectural elements required for sense Alu exonization, identify conserved pre-mRNA structures involved in exon selection and point to a possible chaperone activity of SRP9/14 outside the mammalian signal recognition particle.

Exonic splicing code and protein binding sites for calcium.

Auxilliary splicing sequences in exons, known as enhancers (ESEs) and silencers (ESSs), have been subject to strong selection pressures at the RNA and protein level. The protein component of this splicing code is substantial, recently estimated at ∼50% of the total information within ESEs, but remains poorly understood. The ESE/ESS profiles were previously associated with the Irving-Williams (I-W) stability series for divalent metals, suggesting that the ESE/ESS evolution was shaped by metal binding sites. Here, we have examined splicing activities of exonic sequences that encode protein binding sites for Ca2+, a weak binder in the I-W affinity order. We found that predicted exon inclusion levels for the EF-hand motifs and for Ca2+-binding residues in nonEF-hand proteins were higher than for average exons. For canonical EF-hands, the increase was centred on the EF-hand chelation loop and, in particular, on Ca2+-coordinating residues, with a 1>12>3∼5>9 hierarchy in the 12-codon loop consensus and usage bias at codons 1 and 12. The same hierarchy but a lower increase was observed for noncanonical EF-hands, except for S100 proteins. EF-hand loops preferentially accumulated exon splits in two clusters, one located in their N-terminal halves and the other around codon 12. Using splicing assays and published crosslinking and immunoprecipitation data, we identify candidate trans-acting factors that preferentially bind conserved GA-rich motifs encoding negatively charged amino acids in the loops. Together, these data provide evidence for the high capacity of codons for Ca2+-coordinating residues to be retained in mature transcripts, facilitating their exon-level expansion during eukaryotic evolution.

Development of Multidrug Resistance in Acute Myeloid Leukemia Is Associated with Alterations of the LPHN1/GAL-9/TIM-3 Signaling Pathway.

P-glycoprotein (known as ABCB1 transporter) expression in myeloid blasts of acute myeloid leukemia (AML) or myelodysplastic syndrome (MDS) leads to the commonly observed multidrug resistance. Overexpression of latrophilin-1 was detected in leukemic cells from AML patients. In a previous study, we showed that ABCB1 overexpression is associated with decreased latrophilin-1 expression in MOLM-13/VCR and SKM-1/VCR AML cell variants derived from MOLM-13 and SKM-1 cells by vincristine selection/adaptation. In the present study, we found that if ABCB1 overexpression occurs in myeloid blasts of newly diagnosed MDS patients, latrophilin-1 expression is attenuated. Latrophilin-1 may initiate TIM-3- and galectin-9-mediated immune escape. We demonstrated changes in the expression of both proteins by comparing ABCB1-positive cell variants (MOLM-13/VCR, SKM-1/VCR) with their ABCB1-negative counterparts. Galectin-9 was present in our cell lines in eight protein isoforms for which we identified the respective transcription variants resulting from alternative splicing, and we verified their structure by sequencing. The isoform profile of galectin-9 was different between ABCB1-positive and ABCB1-negative cell variants. The interaction partner of galectin-9 is CD44, and its expression was altered in the ABCB1-positive variants MOLM-13/VCR and SKM-1/VCR compared to their ABCB1-negative counterparts.

Restriction of an intron size en route to endothermy.

Ca2+-insensitive and -sensitive E1 subunits of the 2-oxoglutarate dehydrogenase complex (OGDHC) regulate tissue-specific NADH and ATP supply by mutually exclusive OGDH exons 4a and 4b. Here we show that their splicing is enforced by distant lariat branch points (dBPs) located near the 5' splice site of the intervening intron. dBPs restrict the intron length and prevent transposon insertions, which can introduce or eliminate dBP competitors. The size restriction was imposed by a single dominant dBP in anamniotes that expanded into a conserved constellation of four dBP adenines in amniotes. The amniote clusters exhibit taxon-specific usage of individual dBPs, reflecting accessibility of their extended motifs within a stable RNA hairpin rather than U2 snRNA:dBP base-pairing. The dBP expansion took place in early terrestrial species and was followed by a uridine enrichment of large downstream polypyrimidine tracts in mammals. The dBP-protected megatracts permit reciprocal regulation of exon 4a and 4b by uridine-binding proteins, including TIA-1/TIAR and PUF60, which promote U1 and U2 snRNP recruitment to the 5' splice site and BP, respectively, but do not significantly alter the relative dBP usage. We further show that codons for residues critically contributing to protein binding sites for Ca2+ and other divalent metals confer the exon inclusion order that mirrors the Irving-Williams affinity series, linking the evolution of auxiliary splicing motifs in exons to metallome constraints. Finally, we hypothesize that the dBP-driven selection for Ca2+-dependent ATP provision by E1 facilitated evolution of endothermy by optimizing the aerobic scope in target tissues.

Transposon clusters as substrates for aberrant splice-site activation.

Transposed elements (TEs) have dramatically shaped evolution of the exon-intron structure and significantly contributed to morbidity, but how recent TE invasions into older TEs cooperate in generating new coding sequences is poorly understood. Employing an updated repository of new exon-intron boundaries induced by pathogenic mutations, termed DBASS, here we identify novel TE clusters that facilitated exon selection. To explore the extent to which such TE exons maintain RNA secondary structure of their progenitors, we carried out structural studies with a composite exon that was derived from a long terminal repeat (LTR78) and AluJ and was activated by a C > T mutation optimizing the 5' splice site. Using a combination of SHAPE, DMS and enzymatic probing, we show that the disease-causing mutation disrupted a conserved AluJ stem that evolved from helix 3.3 (or 5b) of 7SL RNA, liberating a primordial GC 5' splice site from the paired conformation for interactions with the spliceosome. The mutation also reduced flexibility of conserved residues in adjacent exon-derived loops of the central Alu hairpin, revealing a cross-talk between traditional and auxilliary splicing motifs that evolved from opposite termini of 7SL RNA and were approximated by Watson-Crick base-pairing already in organisms without spliceosomal introns. We also identify existing Alu exons activated by the same RNA rearrangement. Collectively, these results provide valuable TE exon models for studying formation and kinetics of pre-mRNA building blocks required for splice-site selection and will be useful for fine-tuning auxilliary splicing motifs and exon and intron size constraints that govern aberrant splice-site activation.

Cancer-Associated Substitutions in RNA Recognition Motifs of PUF60 and U2AF65 Reveal Residues Required for Correct Folding and 3' Splice-Site Selection.

U2AF65 (U2AF2) and PUF60 (PUF60) are splicing factors important for recruitment of the U2 small nuclear ribonucleoprotein to lariat branch points and selection of 3' splice sites (3'ss). Both proteins preferentially bind uridine-rich sequences upstream of 3'ss via their RNA recognition motifs (RRMs). Here, we examined 36 RRM substitutions reported in cancer patients to identify variants that alter 3'ss selection, RNA binding and protein properties. Employing PUF60- and U2AF65-dependent 3'ss previously identified by RNA-seq of depleted cells, we found that 43% (10/23) and 15% (2/13) of independent RRM mutations in U2AF65 and PUF60, respectively, conferred splicing defects. At least three RRM mutations increased skipping of internal U2AF2 (~9%, 2/23) or PUF60 (~8%, 1/13) exons, indicating that cancer-associated RRM mutations can have both cis- and trans-acting effects on splicing. We also report residues required for correct folding/stability of each protein and map functional RRM substitutions on to existing high-resolution structures of U2AF65 and PUF60. These results identify new RRM residues critical for 3'ss selection and provide relatively simple tools to detect clonal RRM mutations that enhance the mRNA isoform diversity.

Overexpression of GRP78/BiP in P-Glycoprotein-Positive L1210 Cells is Responsible for Altered Response of Cells to Tunicamycin as a Stressor of the Endoplasmic Reticulum.

P-glycoprotein (P-gp, ABCB1 member of the ABC (ATP-binding cassette) transporter family) localized in leukemia cell plasma membranes is known to reduce cell sensitivity to a large but well-defined group of chemicals known as P-gp substrates. However, we found previously that P-gp-positive sublines of L1210 murine leukemia cells (R and T) but not parental P-gp-negative parental cells (S) are resistant to the endoplasmic reticulum (ER) stressor tunicamycin (an N-glycosylation inhibitor). Here, we elucidated the mechanism of tunicamycin resistance in P-gp-positive cells. We found that tunicamycin at a sublethal concentration of 0.1 µM induced retention of the cells in the G1 phase of the cell cycle only in the P-gp negative variant of L1210 cells. P-gp-positive L1210 cell variants had higher expression of the ER stress chaperone GRP78/BiP compared to that of P-gp-negative cells, in which tunicamycin induced larger upregulation of CHOP (C/EBP homologous protein). Transfection of the sensitive P-gp-negative cells with plasmids containing GRP78/BiP antagonized tunicamycin-induced CHOP expression and reduced tunicamycin-induced arrest of cells in the G1 phase of the cell cycle. Taken together, these data suggest that the resistance of P-gp-positive cells to tunicamycin is due to increased levels of GRP78/BiP, which is overexpressed in both resistant variants of L1210 cells.

Homogentisate 1,2-dioxygenase (HGD) gene variants, their analysis and genotype-phenotype correlations in the largest cohort of patients with AKU.

Alkaptonuria (AKU) is a rare metabolic disorder caused by a deficient enzyme in the tyrosine degradation pathway, homogentisate 1,2-dioxygenase (HGD). In 172 AKU patients from 39 countries, we identified 28 novel variants of the HGD gene, which include three larger genomic deletions within this gene discovered via self-designed multiplex ligation-dependent probe amplification (MLPA) probes. In addition, using a reporter minigene assay, we provide evidence that three of eight tested variants potentially affecting splicing cause exon skipping or cryptic splice-site activation. Extensive bioinformatics analysis of novel missense variants, and of the entire HGD monomer, confirmed mCSM as an effective computational tool for evaluating possible enzyme inactivation mechanisms. For the first time for AKU, a genotype-phenotype correlation study was performed for the three most frequent HGD variants identified in the Suitability Of Nitisinone in Alkaptonuria 2 (SONIA2) study. We found a small but statistically significant difference in urinary homogentisic acid (HGA) excretion, corrected for dietary protein intake, between variants leading to 1% or >30% residual HGD activity. There was, interestingly, no difference in serum levels or absolute urinary excretion of HGA, or clinical symptoms, indicating that protein intake is more important than differences in HGD variants for the amounts of HGA that accumulate in the body of AKU patients.